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We examine the crystallization dynamics of nanoparticles reversibly tethered by DNA hybridization. We show that the crys- 
tallization happens readily only in a narrow temperature "slot," and always proceeds via a two-step process, mediated by a 
highly-connected amorphous intermediate. For lower temperature quenches, the dynamics of unzipping strands in the amor- 
phous state is sufficiently slow that crystallization is kinetically hindered. This accounts for the well-documented difficulty of 
forming crystals in these systems. The strong parallel to the crystallization behavior of proteins and colloids suggests that these 
(3JT) disparate systems crystallize in an apparently universal manner. 
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The use of DNA as a programmable linking agent is a practi- 
cal, "bottom-up" approach to materials design 1 4 . One starts 
from "molecules", consisting of nanoparticles (NP) function- 
alized by multiple single stands of DNA (ssDNA). When the 
DNA on adjacent particles hybridize to form double- stranded 
DNA (dsDNA) the particles are physically linked, potentially 
leading to the formation of complex structures 5 6 . If one can 
specify the number and orientation of the ssDNA attached to 
the NP, it is possible to control the local geometry of the net- 
works, which in turn may control the geometry of higher-order 
structures^. This bottom-up approach can result in a precision 
hard to achieve by molecular nano-fabrication, with promis- 
ing future applications in optical and electrical materials^. 
While there have been some recent successes creating crys- 
talline ordered arrays of DNA-linked NP^^, the formation 
of regularly ordered structures has proved challenging. More 
frequently, the NP assemble into disordered aggregatesEE] 
Therefore, we aim to better understand the dynamical path- 
ways the system must follow in order to crystallize, as well as 
the mechanisms that hinder ordering, with the ultimate goal of 
avoiding kinetic bottlenecks. 

To put these difficulties in the context of more traditional 
materials, experiments on colloidal particles, which isotrop- 
ically interact with each other, are well-known to exhibit a 
two-step crystallization mechanism. In these cases, the gas- 
liquid coexistence curve is metastable relative to the gas- solid 
coexistence curve. Quenching these colloidal system inside 
the gas-liquid coexistence region results in (i) a phase sepa- 
ration into a high-density liquid, followed by (ii) crystal nu- 
cleation within this high-density droplet. Quenching the sys- 
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tern to very low temperature causes it to simply form disor- 
dered gels, kinetically hindering crystal formation. There is 
considerable theoretical and simulation evidence for this sce- 
narioESEl Similar ideas have also been proposed for proteins, 
and it is now believed that there is a crystallization tempera- 
ture "slot" outside which crystallization does not occu r 25 " 28 i 
An alternate mechanism for clustering preceding crystalliza- 
tion has also been proposed for the case of proteins interacting 
through strong, but patchy interactions. It has been argued that 
self-assembly, driven by highly specified local geometry im- 
posed by the bonding sites, can create a locally high-density 
region which enhances the formation of nuclei, even in the 
absence of the thermodynamic drive for phase separation. A 
two-step mechanism, mediated by this self-assembled amor- 
phous state, is thus another pathway by which proteins can 
crystallize 23 . While the factors controlling the crystallization 
of colloids and proteins (or patchy colloids) may be system 
specific, we stress that crystallization in these soft matter sys- 
tems always seems to follow a two-step kinetic scheme, with 
an amorphous, highly-connected phase serving as a kinetic in- 
termediate. 

Motivated by this apparent universality in these systems, we 
examine the crystallization of DNA-linked NP, where the for- 
mation of clusters is controlled by DNA hybridization. We 
show that, in spite of the significant differences in the phys- 
ical connectivity between this case and previously examined 
situations (i.e., colloids and proteins), the crystallization of 
NPs linked by DNA tethers also follows a two-step process: 
initially there is a cluster of linked particles without any crys- 
talline order. This process is facilitated either by self-assembly 
of the nanoparticles, or by phase separation, depending on the 
region of parameter space explored. We show that the persis- 
tence of this amorphous state grows very rapidly on cooling, 
resulting in a very narrow crystallization temperature slot. We 
thus argue that the crystallization of such ssDNA tethered NPs 
follow the same universal behavior as found for colloids, pro- 
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teins and other patchy colloids. 

2 Modeling 

In our simulations, each NP is grafted with 6 chains in an oc- 
tahedral symmetry. This orientation will naturally lead these 
NP to crystallize into a simple cubic structure. We use an ef- 
fective potential model developed to capture the base-pair se- 
lectivity between two ssDNA based on the nucleotides' iden- 
tity (A, T, C, or G), and the bonding specificity which allows 
only one bond to each base 33 34 . The effective potential be- 
tween two ssDNA depends only on the intermolecular sepa- 
ration and their relative angular orientation. The parameters 
of the effective potential are obtained by a systematic coarse- 
graining of a more detailed model for the DNA interactions, 
and it has been verified that the coarse-grained model quantita- 
tively reproduces the behavior of the more explicit model^ES. 
We study this model via Monte Carlo (MC) simulations in the 
canonical ensemble (fixed number of particles N, volume V, 
temperature T). For a randomly chosen NP, we perform three 
sequential attempts - (i) NP translation without rotation, (ii) 
NP rotation without translation, and (iii) a combined NP trans- 
lation and rotation. We define one Monte Carlo step (MCS) as 
N such attempts, where N is the total number of NP. Note that 
the dominant interactions in this system are between the bond- 
ing arms; the short-ranged NP core repulsion plays relatively 
little role. There are no explicit solvent interactions, so we 
cannot directly capture effects such as salt concentration. We 
used reduced units, as defined in ref.^. To improve the statis- 
tics of our results, we average over 5-15 independent runs for 
each system. 

The formation of a crystal or liquid phases in this system 
is driven by the DNA base pairing alone, resulting in very low 
density phases. The addition of significant attractions between 
NP cores - which might be expected in some experimental 
systems - can result in much higher density phases, where the 
density is determined by the packing of the NP cores. Such 
phases are analogous to those found for isotropically interact- 
ing colloids. We do not consider this more complicated situa- 
tion in this paper, and defer it to future research. 

Previously, ref. 36 examined the phase diagram for this 
model, and showed that it exhibits polymorphous phase be- 
havior with at least six distinct crystal phases. The lowest 
density crystal, crystal I, consists of a simple cubic (SC) lattice 
that reflects the octahedral symmetry of the functionalized NP. 
The length of the DNA connections between NP leaves ample 
space, allowing the cubic order to repeat itself as a hierarchy of 
interpenetrating cubic lattices. Accordingly, crystal II consists 
of two interpenetrating SC lattices, crystal III has three inter- 
penetrating SC lattices, etc. Experimentally, achieving this in- 
terpenetration may be challenging, due to repulsions from the 
charged DNA backbone. However, with a proper choice of 
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Fig. 1 The phase diagram of the octahedrally functionalized NP^. 
Solid vertical lines are the densities of crystal I and II, and the im- 
ages show simulation snapshots of the local order of each crystal; 
in these snapshots, the blue spheres represent the core NP units, 
and DNA bonds between NP are represented by the green connec- 
tions. The grey regions of the phase diagram indicate crystal-gas or 
crystal-crystal coexistence. The crystal melts on heating to at the 
bold, dashed line. The faded-red regions are the metastable amor- 
phous phase separation regions. The estimated phase boundaries for 
the gas-liquid (0 < pd 3 < 1) and liquid- liquid (1 < pd 3 < 2) phase 
transitions are indicated by the light dashed lines. The red diamonds 
indicate the state points where we quench to study the crystallization 
dynamics - one point in the phase amorphous separation region (a), 
and two in assembly dominated regimes (b) and (c). These letters 
correspond to those used in Fig. [2] 

solution, electrostatic interactions beyond a few nm between 
DNA can be effectively screened 37 . Accordingly, function- 
alizing strands should be long enough to open the network 
structure and thereby minimize electrostatic repulsions. 

Fig. [T] shows the melting temperature for crystal I and II 
and the metastable amorphous phase boundaries; the phase 
boundary of the metastable amorphous phases is estimated by 
extrapolating the observed phase boundaries of NP functional- 
ized by 3, 4, or 5 DNA strands in refJ^l. Parallel to the crystal 
polymorphism, there is liquid state polyamorphism - that is, 
this system exhibits the unusual feature of multiple thermody- 
namically distinct liquid states in a pure syste m 38 " 40 ! 

We focus primarily on quenches to T — 0.089 where we 
empirically find that crystallization proceeds most readily. We 
also consider other T to test how the rate of crystallization 
varies. We investigate systems at pd 3 = 0.5, pd 3 = 1, and 
pd 3 = 1.5, where pd 3 is the scaled density so that a single SC 
lattice has density pd 3 = 1^1 We keep V fixed, so these three 
densities correspond to N = 500, 1000, 1500, respectively. 
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The system at pd 3 = 0.5 is located in the amorphous phase 
separation region, whereas other two systems at pd 3 = 1, and 
pd 3 = 1.5 are outside of any amorphous-amorphous phase 
transition. These three systems allow us to separately address 
the role of single versus double-interpenetrating crystals, and 
spinodal-assisted versus assembly-driven crystallization pro- 
cesses. 

3 Results 

Since we wish to track the dynamics of crystallization, we 
need to separately identify the formation of clusters and the 
degree of crystallinity of those clusters. First, we evaluate 
the evolution of the cluster size directly from the number of 
bonded NP and calculate the weighted mean cluster size 
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where n is the size of a cluster at some time t, and P(ri) is 
the probability of finding a cluster of size n. Here, a bond is 
determined by the linkage of the DNA strands connecting NP. 

Following the approach in ref. 41 , we can identify crystal- 
like regions using an orientational order parameter g/, defined 
by a sum over spherical harmonic functions !/, w (r) using the 
unit vectors f defined by the bonded neighbors. Among the 
various choices for the degree / of l/ ?m (f), / = 4 provides 
the strongest signal for the expected cubic symmetry; specif- 
ically, Q4 = 0.764 for an ideal SC lattice. To capture the lo- 
cally SC structure in case of two separate lattices, we consider 
only bonded neighbors when calculating Q4, excluding un- 
bonded neighbors, which might have small separation but be- 
long to a distinctly different bonded network. Additionally, we 
compute g6 for nearest neighbors to identify possible body- 
centered-cubic (BCC)-like structures of non-bonded units that 
might arise due to interpenetration at pd 3 = 1, 1.5. For amor- 
phous systems, Q\ = in the thermodynamic limit, so that one 
can immediately distinguish crystal from amorphous systems. 

While Q4 and are useful to identify global crystallinity 
in a given configuration, they are less helpful in identifying the 
presence of local crystal regions embedded in a larger amor- 
phous cluster. Identification of these locally crystalline re- 
gions is necessary to track if crystals can assemble directly, or 
if crystals only form following a connected amorphous inter- 
mediate state. Following ref. 41 , we identify crystal-like par- 
ticles using a local invariant q^{i) for each individual parti- 
cle, and the corresponding complex vector q4 (i) . A particle is 
said to be crystal-like if it has a minimum number of neigh- 
bors with crystal-like connections. Bonded particles i and j 
are said to have a crystal-like connection if the vector dot 
product q4(i) • q4(j) exceeds a threshold value. By compar- 
ing the distribution of dot product values of amorphous sys- 
tems with well-crystallized systems, we find that a threshold 
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Fig. 2 The time evolution of the clustering and crystallization process 
at (a) pd 3 = 0.5, where spinodal decomposition promotes clustering; 
(b) pd 3 = 1, the density of the single SC crystal, where DNA hy- 
bridization drives clustering; (c) pd 3 = 1.5, where there is a phase 
separation of the single and double interpenetrating crystals. For 
each density, we show the fractional cluster size s/N, the fraction of 
crystal-like particles Nx/N, and the normalized average orientation 
G4/G4 C - F° r Pd 3 = 1 and 1.5, we also show Q6/Q^ CC to determine 
whether the formation of two interpenetrating crystals are simultane- 
ous or occur step- wise. 



value q4(i) • q4(j) > 0.95 for the dot product identifies more 
than 90% of bonds in crystal state, and only misidentifies less 
than 2% of the bonds of the amorphous system as crystal- 
like at pd 3 — 0.5. Similar precision is found at pd 3 = 1 and 
pd 3 = 1.5. Finally, analysis of these systems shows that we 
reliably define a crystal-particle if it has at least three crystal- 
like bonded neighbors. We use this as our criterion to label an 
NP as part of a crystal. 

We examine the evolution of the crystallizing systems by 
evaluating the fractional cluster size s/N, the fraction of 
crystal-like nanoparticles Nx/N, and the normalized average 
orientation Qa/Q^ 1 , so that all quantities vary over the range 
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[0, 1] (Fig. [2]). The density pd 3 = 0.5 allows us to examine 
crystallization to a single network in the presence of an amor- 
phous phase separation, similar to the case of colloidal sys- 
tems. For comparison, density pd 3 = 1.0 follows crystal for- 
mation driven only by the assembly of DNA links. Finally, for 
pd 3 = 1.5, self-assembly (without phase separation) drives the 
formation of a higher density state so that the system crystal- 
lizes into a combination of a single cubic lattice and a second 
higher density interpenetrating lattice. Accordingly, we can 
determine if the pathway for crystallization for interpenetrat- 
ing networks differs from that for a single network, and com- 
pare spinodal driven versus self-assembly driven clustering. 

Figure 2 shows that, after quenching from high T to T = 
0.089 (just below the hybridization temperature for the DNA), 
we find that all systems rapidly undergo a condensation from 
an unbonded state to a large, bonded amorphous cluster. Both 
Q4 and Nx remain small for t < 10 4 , demonstrating the amor- 
phous nature of the cluster. The ordering process only occurs 
much later, at t & 10 4 when Q4 and Nx rise sharply within a 
narrow window of time. Since the ordering happens over a 
relatively narrow time window, the global Q4 is itself an in- 
dication of crystal formation. Note that the fraction Nx/N 
never reaches one, since there are always surface particles of 
the crystal that will not be identified as crystal-like. These re- 
sults establish that, for each state point considered, the system 
first forms a highly connected amorphous phase, from which 
a crystal nucleates and grows. We refer to this as the two-step 
process of crystallization. 

Having established the two-step nature of crystallization, 
we next wish to determine if the intermediate amorphous 
phase can be considered a metastable equilibrium, and how 
this impacts the difficulty of the eventual crystallization. To 
do so, we examine the crystallization dynamics as a function 
of the quench depth by evaluating the time Xx needed for the 
crystal to nucleate and the internal relaxation time x a of the 
amorphous intermediate. We define %x by the time when 10% 
of the particles are designated as crystal, since the crystalliza- 
tion process appears irreversible at this fraction; an alternate 
criterion will change the value of Xx, but not the T depen- 
dence. We define x a for the amorphous phase from the re- 
laxation time of the coherent intermediate scattering function 
F(q,t), evaluated at the wave vector q corresponding to the 
bonding distance between NP; we choose this wave- vector 
since it captures the slowest relaxation (apart from q — » 0). 
The determination of x a is complicated due the the aging of 
the amorphous state following the temperature quench. To 
limit aging effects, we wait for the largest possible time to be- 
gin calculating F(q,t) that will still allow F(q,t) to decay to 
zero prior to crystallization. 

Figure [3] shows the time crystallization time Xx and the 
amorphous relaxation time x a for density pd 3 =0.5, which 
have been averaged over 10 independent trajectories. Such 
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Fig. 3 Time-temperature-transformation diagram showing the tem- 
perature dependence of the crystallization time Xx and amorphous 
intermediate relaxation time T a . The minimum of Xx shows there is 
only a narrow slot where crystallization readily proceeds. 

a plot is commonly r eferre d to as a "time-temperature- 
transformation" diagram^E! We find that x a is signifi- 
cantly smaller than Xx so that the clustered state can reach 
a metastable equilibrium prior to crystallization. Hence the 
first step toward crystallization also includes the equilibration 
of the metastable state. On cooling, x a increases rapidly, since 
the lifetime of dsDNA pairs grows quickly, thereby hindering 
relaxation. In contrast, Xx initially decreases on cooling, ex- 
pected since the predicted barrier to crystallization decreases 
on cooling from classical nucleation theory^! However, on 
further cooling, Xx rapidly increases, as it become dominated 
by the slow relaxation of x a . From a physical perspective, 
the persistence of DNA base-pair bonds prevents unzipping 
on a reasonable time scale, and thus the amorphous cluster 
cannot reorganize to "find" the crystal state, resulting in a ki- 
netically dominated process. Accordingly, there is a very nar- 
row slot that must be found for successfully nucleating the 
crystal state. The minimum, or "nose", in the crystallization 
time is ubiquitous in supercooled liquids^^. If the system 
is cooled below r nose in a time less than T n0S e> the system will 
not have adequate time to crystallize. Accordingly, the nose 
defines a critical cooling rate R= (Tm — T n0SQ )/x n0S Q, where 
Tm is the melting temperature 42 ; cooling faster than this rate 
will prevent crystallization in all cases, and should therefore 
be avoided. 

For the highest density pd 3 = 1.5, we wish to further deter- 
mine if the formation of a double network occurs simultane- 
ously with, or after the formation of a single SC network. To 
test this, we check for the presence of BCC order in the lattice 
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using 26 (since it is more sensitive to BCC order) and compare 
its evolution to Q4 . For the calculation of g6 , we use spatial 
separation, rather than bonds, to determine neighbors, since 
the units comprising the BCC structure are actually unbonded 
neighbors in a separate cubic lattice. These unbonded neigh- 
bors are separated by a distance less than the bonding dis- 
tance. Choosing a cutoff separation of V?>/2d w 0.866 d (ra- 
tio of BCC to FCC lattice spacing) effectively excludes cubic 
bonded neighbors, and includes most neighbors that should 
have BCC order. We actually use a slightly less restrictive 
definition, and include neighbors up to a distance 0.92 d, since 
the positions of interpenetrating particles are not rigidly fixed; 
this cutoff is still small enough that it excludes the vast major- 
ity of bonded neighbors. 

For reference, we first examine density pd 3 = 1.0 where 
there should be no interpenetration, so that we know to 
what degree (?6 might give a false signal of interpenetration. 
Fig^b) shows that Q 6 /Q^ CC ~ 0.3 for the final SC crystal, so 
that we have little false signal of interpenetration. The small 
value of Qe / Q^ cc can be largely attributed to our choice of 
the cutoff distance used to determine neighbors, since if we 
included the bonded neighbors of the SC lattice we would 
expect Q 6 /Q% cc = 0.693 (for a perfect SC lattice). Apply- 
ing this metric for the case with interpenetration (pd 3 = 1.5), 
Fig.^c) shows that Qe/^ CC captures the eventual interpen- 
etration. However, the growth of g6 slightly lags behind that 
of of Q4. Thus, the formation of a single cubic lattice appears 
to slightly precede the interpenetrating structure. Note that the 
asymptotic value of Nx/N for pd 3 = 1.5 is slightly larger than 
for the lower densities. This is because we use a larger N, and 
thus the surface-to- volume ratio is smaller, and so the surface 
effect on Nx/N is smaller. 

For both densities, note that s slightly increases upon crys- 
tallization. Presumably, the long ranged order allows for the 
formation of some additional bonds that were dangling in the 
clustered, but disordered state. For pd 3 = 1.5, s has a very 
small but noticeable decrease immediately prior to ordering. 
This can be understood by the fact that the formation of two in- 
terpenetrating SC lattices requires a separation of two lattices; 
therefore, there must be a temporary breaking of amorphous 
bonds between locally cubic lattices before ordering, which 
temporarily decreases s. 

We directly visualize the crystallization of two systems 
(pd 3 = 0.5 and pd 3 = 1.5) at three important points: (i) clus- 
tering prior to ordering; (ii) the ordering process, and (iii) after 
ordering (Fig. [4]). The crystal-like NP are colored in red with 
a slightly larger radius. From Fig.|4](b) and (e) we can clearly 
see the interface between crystal and amorphous phases. Fig. [4] 
(b) has two separate SC lattices and evidence of multiple nu- 
cleation cores. 
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Fig. 4 Visualization of the system at three points during crystalliza- 
tion. (a)-(c) are at density pd 3 =0.5. Specifically, (a) is system after 
clustering but prior to ordering; (b) is during the ordering; (c) is after 
ordering, (d)-(f) represent the same time progression but at density 
pd 3 = 1.5. The crystal-like particles are colored in red with a slightly 
larger size. 



4 Discussion and Conclusions 

The observed sequence of clustering via DNA links followed 
by ordering for all densities confirms that the crystalliza- 
tion dynamics for the DNA-linked NP follows a conventional 
"two-step" pathway of crystallization - even in the case where 
a double-interpenetrating network must form. Depending on 
density, the amorphous intermediate is driven either through 
metastable phase separation or DNA driven assembly. Addi- 
tionally, the intermediate amorphous state has a rapidly grow- 
ing lifetime on cooling. Thus, the same framework used to 
understand protein and colloid crystallization kinetics also ap- 
plies to this more unusual material. 

We compare our results with the recent experimental stud- 
ies of Mirkin and co-workers 14 , where they report a 3 -stage 
crystallization process for uniformly DNA-coated NP. In that 
study, the NP initially form small amorphous aggregates, and 
these clusters separately evolve crystallinity at stage two. The 
ordered but dispersed clusters eventually coalesce into a large 
final crystal lattice. As noted by these authors, this last stage 
of crystallization is probably driven by Ostwald ripening, a 
mechanism that is relatively well explored. Accordingly, these 
experiments also fit within the general framework of the two- 
step process, since the crystallinity of small clusters evolves 
from the amorphous aggregates. To further compare to these 
experiments, we also checked if very slow ripening might oc- 
cur in our simulations under different thermodynamic con- 
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ditions. Indeed, simulations at lower density (pd 3 = 0.3) 
evolve crystallinity in small region within the amorphous clus- 
ter. The subsequent growth of this crystal is so slow, that we 
are not able to complete the crystallization in the computa- 
tional time frame of the simulation. This slower scenario for 
crystal growth is likely the same as that observed in ref H Ac- 
cordingly, the rate of the "second" step in the process (crystal 
growth) can vary significantly depending on the state point 
chosen. We should note a potentially important difference be- 
tween our simulations and that of the experiments of ref. 14 . 
For the uniformly coated NP used experimentally, the local 
NP binding is not predisposed to reflect the order of the even- 
tual crystal state. In contrast, our 6-armed units bind in such a 
way to readily reflect SC symmetry. We expect that the coars- 
ening to the eventual crystal state should occur more readily 
in the case where the crystal order is reflected in the symmetry 
of local bonding. Nonetheless, since the experimental system 
and the simulated model behave similarly, we conclude that 
the two-step crystallization is the apparently universal path- 
way followed by these systems. 

Experiments using nanoparticles functionalized with many 
strands of DNA (as opposed to a small number of strands) 
present an additional potential barrier to crystallization. As 
the strand density increases, it has been experimentally found 
that the hybridization transition becomes increasing ly sharpie. 
This behavior has also been argued for theoretically^! Nar- 
rowing the hybridization window will also result in narrowing 
the crystallization temperature slot - making the formation of 
ordered structures even more challenging. Accordingly, it may 
be valuable to work with NP functionalized by a small number 
of strands. These limited functionality systems also have the 
advantage that they lend themselves more readily to a theoret- 
ical description 47 . 

Since there are no NP attractions in our model (other than 
indirectly via DNA linking), the gas-liquid phase separation 
is driven by the self-assembling DNA-hybridization. If, in 
addition to the hybridization, there were isotropic NP attrac- 
tions, the system might be able to form a much higher density 
droplet controlled by the packing of the NP cores. This could 
result in a crystallization by hybridization that would actually 
dramatically decrease the density relative to the amorphous 
cluster, since DNA links will serve to open the structure, qual- 
itative similar to the crystallization of ice from water. Such an 
intermediate not created by DNA base pairing might offer dif- 
ferent pathways to creating crystals, without the kinetic traps 
that are encountered experimentally. This is another possible 
avenue to pursue to experimentally facilitate the crystalliza- 
tion process. 

The highly specified DNA orientation in our model clearly 
helps to build in self-assembly of the higher-order crystal 
structure. However, this is not a sufficient condition to guar- 
antee the ready formation of crystals. If that were the case, a 



previous study of NP decorated with ssDNA in a tetrahedral 
orientation 35 should have discovered the spontaneous forma- 
tion of diamond lattices, but instead found only amorphous 
gels. Hence, significant work remains to understand how to 
best design a desired higher order structure from relatively 
simple building blocks. 
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